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Abstract 



We examine Gamow's method for calculating the decay rate of a wave 
function initially located within a potential well. Using elementary techniques, 
we examine a very simple, exactly solvable model, in order to show why it is 
so reliable for calculating decay rates, in spite of its conceptual problems. We 
also discuss the regime of validity of the exponential decay law. 
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I. INTRODUCTION 



Complex-energy "eigenf unctions" made their debut in Quantum Mechanics through the 
hands of Gamow, in the theory of alpha-decay [|l|. Gamow imposed an "outgoing wave 
boundary condition" on the solutions of the Schrodinger equation for an alpha-particle 
trapped in the nucleus. Since there is only an outgoing flux of alpha-particles, the wavefunc- 
tion ipi^yt) must behave far from the nucleus as (for simplicity, we consider an s-wave, and 
use units such that h = 2m = 1) 

^—lEt+ikr 

ipij-.t) (r oo). (1) 

This boundary condition, together with the requirement of finiteness of the wave function 
at the origin, gives rise to a quantization condition on the values of k (and, therefore, on 
the values oi E = k'^). It turns out that such values are complex: 

kn = l^n ~ "^-^n/S, E^ = — i^n/'^) (2) 

and so it follows that 

-r.at+Kr,r 

(r-.oo). (3) 

Thus, if Vn > 0, the probability of finding the alpha-particle in the nucleus decays exponen- 
tially in time. The lifetime of the nucleus is then given by r„ = l/r„, and the energy of the 
emitted alpha-particle by e„. 

Although very natural, this interpretation suffers from some difficulties. How can the 
energy, which is an observable quantity, be complex? (In other words, how can the Hamil- 
tonian, which is a Hermitean operator, have complex eigenvalues?) Also, these "eigenfunc- 
tions" are not normalizable, since r„ positive implies Kn positive and so, according to (^, 
\ipn{r,t)\'^ diverges exponentially with r. 

In spite of such problems (which, in fact, are closely related), it is a fact of life that 
alpha-decay, as well as other types of decay, does obey an exponential decay law and, in 
many cases, Gamow's method provides a very good estimate for the decay rate. Why this 
method works is a question that has been addressed in the literature using a variety of 
techniques ■ Here we examine this question in a very elementary way, using techniques 
that can be found in any standard quantum mechanics textbook and some rudiments of 
complex analysis. 

This paper is structured as follows. In Section we show Gamow's method in action for 
a very simple potential. Some of the results obtained there are used in Section |Tl|, where we 
study the time evolution of a wave packet initially confined in the potential well defined in 
Section ||. This is done with the help of the propagator, built with the true eigenf unctions 
(i.e., associated to real eigenenergies) of the Hamiltonian. As a bonus, we show that the 
exponential decay law is not valid either for very small or for very large times. This is the 
content of Section where the region of validity of the exponential decay law is roughly 
delimited. 
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II. DECAYING STATES 



In order to exhibit Gamow's method in action, we shall study the escape of a particle 
from the potential well given by: 

Y(^) = I -a) fora;>0, .^s 

^ ^ I +00 for X < 0. ^ ^ 

(Escape from this potential well was studied in detail in Refs. [^,^,10|. In this section we 
follow closely the treatment of Ref. .) The positive dimensionless constant A is a measure of 
the "opacity" of the barrier at x = a; in the limit A 00, the barrier becomes impenetrable, 
and the energy levels inside the well are quantized. If A is finite, but large, a particle is not 
confined to the well anymore, but it usually stays there for a long time before it escapes. If 
A is small, the particle can easily tunnel through the barrier, and quickly escape from the 
potential well. Metastability, therefore, can only be achieved if the barrier is very opaque, 
i.e., A 1. Since this is the most interesting situation, we shall assume this to be the case 
in what follows. 

To find out how fast the particle escapes from the potential well, we must solve the 
Schrodinger equation 

d X 
i — i:{x, t) = -— ij{x, t) + -6{x-a) ij{x, t). (5) 

ip{x,t) = exp{—iEt) ip{x) is a particular solution of this equation, provided ip{x) satisfies 
the time-independent Schrodinger equation 

(P A 

- -^ip{x) + -6{x -a)(p{x) = Eip{x). (6) 

Denoting the regions < x < a and x > a by the indices 1 and 2, respectively, the 
corresponding wave functions (pj{x) (j = 1, 2) satisfy the free-particle Schrodinger equation: 

Since the wall at the origin is impenetrable, ipi{0) must be zero; the solution of Eq. (|^ 
which obeys this boundary condition is 



ipi{x) = A sinkx {k = vE). (8) 

To determine (p2{x), we follow Gamow's reasoning [|I|,|^,[TT]| and require ^P2{x) to be an out- 
going wave. Therefore, we select, from the admissible solutions of Eq. (^, 

Mx) = Be'''^. (9) 

The wave function must be continuous at x = a, so that (pi{a) = ^2(0) , or 

— = e-''"' sin ka. (10) 
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On the other hand, the derivative of the wave function has a discontinuity at x = a, which 
can be determined by integrating both sides of (P) from a — e to a + e, with e — ^ O"*": 

- ^[(a) = - (p2ia), (11) 

from which there follows another relation between A and B: 

B kae~^^'^ cos ka 



(12) 



A X — ika 

Combining ([lOD and (|T^) , we obtain a quantization condition for k: 

ka cotan ka = —X + ika. (13) 



The roots of Eq. (|T3D are complex and situated in the half-plane Im/c < 0; when A ^ 1, 
those which are closest to the origin are given by |^,^ 

, mrX ■ fn-nX"^ , ^ ^ , , ,^ 

^^"^ - YTx~ \~X) (^ = l'2,---;^7r<A). (14) 

(For each one of these roots, which are located in the fourth quadrant of the complex A;-plane, 
there is a corresponding one in the third quadrant, given by — fc*. The latter are associated 
to "growing states" and play no role in what follows.) The corresponding eigenenergies 
are 

(15) 

The imaginary part of En gives rise to an exponential decay of |^n(a;, t)P, with lifetime equal 

to 

rn = l/r„ ^ (16) 

Since the corresponding value of B/A is very small (~ n/X), one may be tempted to say 
that the probability of finding the particle outside the well is negligible in comparison with 
the probability of finding the particle inside the well. Normalizing ip^ in such a way that 
the latter equals one when t = 0, the probability of finding the particle inside the well at 
time t, if it were in the n-th decaying state at t = 0, would be 

P„(t)= r\Mx,t)\^dx = exp{-Tnt). (17) 
Jo 

The trouble with this interpretation is that Im A;„ = —Kn/2 < 0, and so ipn{x, t) diverges 
exponentially as x — oo, since, according to (P), 

|V'„(x,t)P = exp(-r„t + K„x) (18) 

outside the well. Because of this "exponential catastrophe", the decaying states are non- 
normalizible and, therefore, cannot be accepted as legitimate solutions of the Schrodinger 
equation (although one can find in the literature ||12[ the assertion that they are "rigorous" 
solutions of the time-dependent Schrodinger equation). 
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III. TIME EVOLUTION OF A WAVE PACKET 



We now return to Eq. (|^) and write, for the solution in region 2, instead of the sum 
of an outgoing plus an incoming wave: 

(^2(x) =e-*'=^' + 5e*'=^ (19) 

Continuity of the wave function at a; = a implies 

A sin ka = e-''"' + Be'^"". (20) 

As before, the derivative of the wave function has a discontinuity at a; = a, given by Eq. 
(p!T|), from which it follows, instead of (f2|). 



kA cos ka = -{- + ik] e-'''^ - ( - - zA; ) 5 e''"'. (21) 



a 



Solving (gOD and (|21|) for A and 5, we find 



2ika 

ka + X e^^"- sin ka ' 
ka + X e"''"' sin ka 



m = -IT^^TI^^r— (22a) 



B{k) = - , ' ■ , ■ (22b) 

^ ^ ka + X e'-'"' sm ka 

These expressions show a couple of interesting features: 

• \B\ = 1 for real values of k, implying a zero net flux of probability through x = 
a; therefore, unlike the solution found in the previous section, there is no loss or 
accumulation of probability in the well region. 

• |A| <^ 1 if ka <^ A, except if k is close to a pole of A{k), in which case \A\ may become 
very large. 

To find the poles of A we must solve the equation A{k)^^ = 0, which, after some algebraic 
manipulations, reads 

ka cota.n ka = —X + ika. (23) 

This is the same as Eq. ([T3|) ! Is this a coincidence? In fact, no. According to (p2|), A and B 
have the same poles; therefore, near a pole both \A\ and \B\ are very large, and Eqs. ( pO[ ) 
and (pID become equivalent to Eqs. ([I0|) and (|T2D , respectively. In what follows, we shall 
show that the poles of A (and B) play an important role in the decay process. 

Suppose that at t = the particle is known to be in the region x < a with probability 
1; in other words, its wave function tp{x, 0) is zero outside the well. Then, at a later time t, 
the wave function is given by 

^{x, t) = r G{x, x'; t) ij{x', 0) dx', (24) 
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where the propagator, G{x,x';t), can be written as 

poo „ 

G{x, x'- 1) = / e-''' * ^fc(x) ifilix') dk. (25) 
Jo 

The function (pki^) is the solution of Eq. (H) corresponding to the energy E = k"^: 

/ ^_ 1 ( A{k) sin kx for x < a, , . 

'^^■^'^^ ~ 71^ ^ I e-*'^" + 5 (A;) e^^^" for x > a. ^ ^ 



With this normahzation, these functions satisfy the completeness relation |13| 

x> 

y^,{x)y,l{x')dk = 6{x-x'). (27) 
Eqs. (|2|)-(gD give, for x < a, 

^(x,t) = — / e-'^ * 0(fc) I A(fc) |2 sin kx dk, (28) 

ZTT JO 

where 

(f){k) = r ^{x',0) sin kx' dx'. (29) 

It is clear that the integral over k is dominated by the resonances, i.e., the neighborhood of 
the poles of A{k). 

Since, for t > 0, e~*^ * —>■ when \k\ ^ oo in the fourth quadrant, one can rotateQ the 
integration contour by 45° in the clockwise sense (see Fig. 1), thus obtaining 

/* OO ^^^^ 

^{x, t) = e"*"/^ / e"^'* /(e"''^/^ k, x)dk+Y\ C{kn, x) e-'^-\ (30) 

•^0 n=l 



where 



and 



f{k,x) = ^(l){k)\A{k)\^ sin kx (31) 



C(fc„, x) = -2m lim {k - fc„) f{k, x). (32) 

The sum in ( |5(]| ) takes into account the poles of A{k) which are situated in the region 
— 7r/4 < argA; < 0, and it corresponds to an expansion in Gamow states (for x < a). 

Let us put aside, for a moment, the integral in (^) (it will be discussed in the next 
section). Then, the "nonescape" probability (i.e., the probability of finding the particle 
inside the well) is given by 



^For this to be possible 4'{k) must be an analytic function of k, but this can be shown to be the 



case |14| if V'(ic',0) is continuous in [0, a]. 
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P{t) = I \il){x, t)p (ix ~ Cn e"'""* + interference terms, (33) 
-'^ n=l 

where c„ = /q" |C(A;„, x)p (ix. For A 1, the interference terms are usually negligible, for 
kn ~ rivr/a and, therefore, the functions C{kn,x) oc sinA;„x are approximately orthogonal. 
On the other hand, since the decay rate r„ of the n-th decaying mode is a rapidly increasing 
function of n (Tn ~ Fi), the decay becomes almost a pure exponential one when Fit > 1. 
The system, therefore, "loses memory" of the initial state. 

Finally, let us note that no exponential catastrophe occurs with ip{x,t). In fact, one can 
easily show, using (0), (^), (0) and the orthogonality of the eigenf unctions v^fc(a;), that 

roo poo 

/ Wx.t)? dx = Wx,Q)?dx, (34) 
JO Jo 

so that an exponential growth of \ip{x,t)\^ is completely ruled out. 



IV. BREAKDOWN OF EXPONENTIAL DECAY 



In order to derive expression (|33D for the nonescape probability, we had to neglect the 
first term on the right hand side of (pO]). In this section we show that such approximation 
is not valid either for very small or for very large times. That it cannot be valid for very 
small t follows from the fact that initially the decay is slower than exponential H,H]. This 



can be easily proved with the help of the continuity equation [|T5|] : 



d , ^ h 

— Pit) = Im 

dt ^ ^ m 



d 



(35) 



Since, by hypothesis, ip{a,0) = 0, it follows that dP/dt = when t = 0, whereas for the 
expression (|33D one has dP/dt ~ — I] CnF^ < at t = 0. 

On the other hand, the exponential decay does not last forever. After a sufficiently long 
time, it obeys a power law P-|5|,|8|,p!6[] . To see this, note that the integral in (|30|), which we 
shall denote here by I{x,t), is dominated by small values of k when t ^ oo, and so can it 
be approximated by 

I{x,t)^^<P'iO)\Am'x k'e-'^'dk^^ (36) 

Therefore, the nonescape probability behaves asymptotically as0 

P(t)^£|/(x,t)pdx~^. (37) 



^Garcia-Calderon, Mateos and Moshinsky 
as when t — > oo, in contrast to Eq. (| 
properly corrected |17|, it also leads to P{t) ~ asymptotically. 



argue that the nonescape probability P(t) decays 
However, there is a flaw in their argument; when 



7 



Comparing (|37D with (0), and using ([I6D, one finds that they become comparable in mag- 
nitude when 

or, since A ^ 1, when 

— ~ 10 In A. (39) 

n 

In practice, when the decay begins to obey a power law the nonescape probability is so small 
(~ A~^°) that it should be very difficult to observe deviations from exponential decay. (On 
the other hand, experimental evidence for non-exponential decay at small times has been 
reported recently [|T8|] .) 

In closing this section, we would like to remark that the breakdown of the exponential 
decay law for either small or large times is not a peculiar feature of the potential (§). It is 
possible to show that an exponential decay cannot last forever if the Hamiltonian is bounded 
below [pJj^^, and cannot occur for very small times if, besides that, the energy expectation 
value of the initial state is finite — conditions which certainly must be satisfied by any 
realistic Hamiltonian or state. 

V. CONCLUSION 

In this paper we showed that decaying states, although plagued by the exponential 
catastrophe, give a fairly good description of the decay of a metastable state, provided 
some conditions are satisfied. In fact, the main objective of this paper was to show that 
one can compute the decay rate solving the time independent Schrodinger equation subject 
to the "outgoing wave boundary condition." This is far from being a trivial result, since 
the corresponding eigenstates are unphysical. The effectiveness of the decaying states in 
describing the decay may be understood by noticing |^ that they are good approximate 
solutions to the time-dependent Schrodinger equation, although nonuniform ones (i.e., they 
are not valid in the entire range of values of t and x). 
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FIGURES 
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FIG. 1. Complex fc-plane. The poles of A{k) are represented by the small circles. Those 
in the fourth quadrant give rise to the sum over decaying modes in Eq. (|3^) when one rotates 
the integration contour of Eq. ( pS]) — the positive real semi-axis — by 45° in the clockwise 
sense (dashed line). 
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